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ABSTRACT 

Although General Relativity (GR) has been tested extensively in the weak grav- 
ity regime, similar tests in the strong gravity regime are still missing. In this paper 
we explore the possibility to use X-ray spectropolarimetric observations of black holes 
in X-ray binaries to distinguish between the Kerr metric and the phenomenological 
metrics introduced by Johannsen and Psaltis (2011) (which are not vacuum solutions 
of Einstein's equation) and thus to test the no-hair theorem of GR. To this end, we 
have developed a numerical code that calculates the radial brightness profiles of accre- 
tion disks and parallel transports the wave vector and polarization vector of photons 
through the Kerr and non-GR spacetimes. We used the code to predict the observa- 
tional appearance of GR and non-GR accreting black hole systems. We find that the 
predicted energy spectra and energy dependent polarization degree and polarization 
direction do depend strongly on the underlying spacetime. However, for large regions 
of the parameter space, the GR and non-GR metrics lead to very similar observational 
signatures, making it difficult to observationally distinguish between the two types of 
models. 

Subject headings: accretion, accretion disks, black hole physics, gravitation, polariza- 
tion, relativity, X-rays: binaries, instrumentation: polarimeters 



1. Introduction 



The theory of General Relativity (GR) has been tested with high accuracy in the regime of 
weak gravity including tests within our solar system and tests based on observations of binary 
pulsars and the double pulsar (|Willlll993l . 120061 ). Although X-ray energy spectra of mass accreting 
stellar mass and supermassive black holes probe GR in the strong-gravity regime, the statistical and 
systematic uncertainties are sufficiently large that the data colle cted so far constrain the parame- 
ters describing the observed systems (i.e. the black hole spin, s ee lMillerll2007l ; iRoss Fabianll2007l ; 



McClintock et al. 



2011 



Kulkarni et al. 



2011 



Gou et al. 



2011 



but do not yet permit to perform 
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sensitive tests of GR (however, see lBambi Baraussdl201ll ; |Bambj |2012j) . In the next few years, 
it may become possible to confront GR with new types of experimental data. In this paper, we 
discuss non-imaging spectropolarimetric X-ray observations of s tellar mass black h oles enabled by a 
mission like GEMS (Gravity and E xtreme Magnetism SME X) (jBlack et alJboid ) or BEST (Black 
Hole Evolution and Space Time) (jKrawczvnski et al.ll2012f ). Other electromag netic observations 
with the potential to test strong gravity GR include pul sar timing observation s (jWex Kopeikin 
1999), radio imaging of supermassive black holes (e.g. iDoeleman et al.ll2009i ). time resolved X- 
ray observations of the Fe K-q fluoresce nt line from galactic and extragalactic black holes (e.g. 
Reynolds &: Nowa k 2003; iGuainazzil |2009| ) , and the observation of stars orbiting the supermassive 
black hole at the center of the Milky Way (jWilll 12008 ; iMerritt et al.ll2010l ). The continued improve- 
ment of the sensitivity of ground-based gravitational wave detectors should enable us to observe 
gravitational waves from compact object mergers before the end of this decade (e.g. lHughes!l2010l . 
and references therein). The ring down following a merger event contains information about dy- 
namical GR in the strong gravity regime, but high signal to noise ratio observations are needed to 
measure the ring down with the required accuracy. 

Observations of phenomena close to the event horizon of stellar mass black holes can probe 
gravity in more than five order of magnitude deeper potential wells, and for scal ar curvatures 



more than twelve orders o f magnitude larger than the best observations today (see iPsaltis 



2008 



Psaltis &; Johannsen 



2011 



, for recent reviews of strong GR tests). GR makes clear predictions of 
how astrophysical black holes look like: they are described by the Kerr metric, a family of vacuum 
solutions that depends on two parameters, the mass M of the black hole, and the angular mo- 
mentum per unit mass a with < a < 1 (in geometric units) corresponding to astrophysical black 
holes, while a > 1 gives rise to a naked singularity and is believed not to correspond to physical 
realizations. A series of papers in the late sixties and seventies established the "no-hair-theorem" 
stating that the Kerr (and Kerr-Newman) family of solutions are th e only stationary axially sym- 
metric vacuum solutions of Einstein's equation (see iRobinsonl |2009| . and references therein). One 
way of testing strong gravity GR thus consists in testing if the Kerr solution indeed describes as- 
trophysical black holes. In this paper, the possibility to use spectropolarimetric observations of 
the X-ray emission from black holes in X-ray binaries to verify GR in the strong gravity regime 
is discussed. For moderate accretion rates, the disks are believed to be thin and to be described 



to good approximation by the Novikov-Thorn e equations (INovikoy &; Thorndll973l ). the relativis 



tic version of the Shakura-Sunyaev equations (jShakura feSunyaevlll973l ). True disks are probably 



somewhat brighter close to the innermos t stable circular orbit (ISCO) of the a ccreting matter owing 



to a non-vanishing torque at the ISCO (jNoble et al. 



2011 



Penna et al.1 12012 ). The radiation from 



a black hole in the "thermal state" is strongly dominated by the thermal em ission from the inner 

and presents a good 



accretion disk (Remillard fc McClintock 



2006 



Abramowicz &; Fragile 



2011 



opportunity to test the validity of accretion disk models and possibly also to test the underlying 
spacetime. 



Quantitative tests of GR should give limits on parameters quantifying the deviation of reality 
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from the predictions of GR. In the weak-gravity regime, a large body of experimental data has been 
used to constrain the deviations of the parameters of the parametrized post-Newtonian formalism 
from their values predicted by GR. In the absence of a general parametrization of the deviation 
from the Kerr spacetime, one has to turn to specific alternatives. Bl ack hole metrics following from 
alternative theories of gravity have been discussed in the literature. lAliev k Gumriikcuoglul (|2005) 
explored a metric describing charged spinning black holes localized on a 3-brane in the Randall- 
Sundr um braneworld. Slowly and rapidly spinning black holes i n Che rn-Si mon gravity were stud - 
ied bvlYuries k Pretoriusl (fcood ) . iKonno. Matsuvama k Tandal (fcood ). and kleihaus et al.1 (|201ll ). 
Pani et al.l (|201ll ) considered slowly rotating black holes in theories where the Einstein-Hilbert ac- 
tion is supplemented by a scalar field coupling to the quadratic algebraic curvature invariants. In 
this paper, we follow the pragmatic approach to compare the p redictions from the K e rr me tric with 
the predictions from the axially symmetric metric described bv lJohannsen k Psaltisl (|201ll ) (JP11). 
The metric includes the Kerr metric as a limiting case, and depends on parameters that describe 
the deviation from the Kerr metric. It can describe static and rapidly spinning black holes and 
does not exhibit pathologies like timelike closed loops outside the event horizon. The objective of 
this paper is to see if X-ray polarimetry has sufficient diagnostic power to constrain the deviations 
from the Kerr metric. 

A number of authors have studied the polarization properties of the X-ray emission from accret- 
ing black holes. The X-rays from a fl at-space Newtonian a c cretion disk are expected to be polarized 
owing to Thomson scatter in g (e.g. 



Chandrasekhar 



1960 



Angell Il969l: lLoskutov k Sobolev 



1982: 



Sunvaev k Titarchukll 19851 ). IStark k Connors! (|1977l ). IConnors et al.l (119801 ). and lConnors k Stark 
(1980) showed that accounting for GR effects, the polarization degree and polarization direction of 



the emission from a thin accretion disk show a complex energy dependence which might be used to 
estimate the pa rameters descr i bing t he system, e.g. the inclination of the system and the mass of 
the black hole. IDovciak et al.1 (120081) evalua t ed th e impact of various atmospheric optical depths 
on the observed polarization signal. iLi et al.l (|2009l ) studied the same problem with a focus on how 
the polarization information can b e used to break model degen eraci es, i. e. to constrain th e incli - 



Matt. Fabian, k Ross 



Poutanen k Svensson 



(|1993h 
199&h 



nation of the inner accretion disk. lLaor. Netzer. k Piranl (|l990l ) and 
analyzed the polarization of the UV/soft X-ray emission from AGNs. 
calculate d the polarization of em issi on which Comp t on scatters in a hot corona above the accretion 

) discussed the polariza tion of emis- 



disk, and lPoutanen et al.l (j 19961 ) and 



sion scattered by accretion disks. IDovciak et al 



Dovciak et al 



(2004, 
(l2Q0jand 



2011 



Zamaninasab et a 



(2011 



scrutinized 



the po lariz ation signature produ ced by orbiting hot-spots. iDavis et al.1 (|2009l ). ISilant'ev k Gnedin 
(120081 ) and ISilant'ev et al.1 (120111) studied the impac t of turbulent magnetic fields on the observed 
polarization signatures. ISchnittman k Krolikl (|2009l ) emphasized the importance of radiation scat- 
tered after returning to the disk owing to the spacetime curvature in the surrou nding of the black 
hole on the observed polarization (see ICunninghamlll976l ; lAgol k Krolikl |2000| . for earlier related 
studies). The same authors studied the observational app earance of the hard and ste ep power-law 
(SPL) states when coronal emission cannot be neglected (jSchnittman k Krolikll2010l ). 
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In this paper, we extend the work of lSchnittman Krolikl (|2010l ) to cover not only Kerr space- 
times but also the metric of JP11. Section [2] will discuss the technical aspects of our calculations, 
including the equations used to derive the radial structure of the accretion disk, the formalism 
adopted to calculate the polarization of the emission, and the ray tracing code. Section [3] will de- 
scribe the results of the simulations including the predicted observational signatures of the non-GR 
metrics. We end with a discussion of the results in Sect. [U 



We use the notation of iMisner. Thorne. fc Wheeled (|1973l . MTW73). Distances are given in 



units of the gravitational radius r g = G M/c 2 , and we set G = c = K = 1 throughout the paper. 
Einstein's summing convention is used. Greek indices run from to 3 and Latin indices from 1 to 
3. 



2. Methodology 

2.1. The Metric, Stable Circular Orbits, and Transformation Matrices 

The following treatment assumes that Einstein's Equivalence Principle (EEP) is valid. In 
brief: the trajectories of freely falling "test" bodies are independent of their internal structure 
and composition; the outcome of any local non-gravitational experiment is independent of the 
velocity of the freely-falling reference frame in which it is performed; the outcome of any local 
non- gravitation al experiment is independent of where and when in the Universe it is performed 



(see IWhlll2006l ). We adopt the postulates of metric theories (the only theories which are consistent 
with the EEP) namely that (i) spacetime is described by a symmetric metric, (ii) test bodies follow 
geodesies of this metric, and (iii) the non- gravitational laws of physics are consistent with the theory 
of special relativity in local freely falling reference frames. 

The metric of JP11 reads: 

,2 r, w «m A 2Mr\ , 2 4aMrsin 2 H ,. , „ E[l + h(r,9)] r2 
ds 2 = - [1 + h(r, 0)] 1 - -=- dt 2 [1 + h(r, 6)]dtdtf> + - 9 . 2 ' / ~ dr 2 



+ Ed9 2 + 



2 f 2 2 2a 2 Mrsin 2 9\ o 2 (E + 2Mr) sin 4 9 



shr 9 \ r z + a z H + h(r, 



(1) 



with E = r 2 + a 2 cos 2 9 and A = r 2 — 2Mr + a 2 . The constants M and a are the mass and the spin 
of the black hole as measured by an observer at infinity. 

The function 



h(r,9)=e 3 ^ 2 - (2) 



M 3 r 

parametrizes the deviation from the Kerr metric to leading order, as lower orders are tightly con- 
strained by the requirement for asymptotic flatness, the absence of gravitational waves, and Lunar 
Laser Ranging results (see JP11). For 63 = the metric reduces to the Kerr metric in Boyer 
Lindquist coordinates. We call = (t, r, 9, <j>) the Global Coordinates (GCs) in the following, and 
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denote the basis vectors of the tangent vector space at a point V as = d/dx^. The largest root 
of the equation g t ? — gug^ = gives the location of the event horizon, and we infer the location 
of the ISCO from the condition: dE/dr = 0. 

The Killing vectors and lead to the conserved energy E = —pt and angular momentum 
L z = at infinity. JP11 derived these two quantities for particles of mass /x orbiting the black hole 
on circular geodesic orbits at radial coordinate r. The equations read (without loss of generality, 
we set {j, = 1): 

1 /Pi + P 2 

(3) 



Pa 



L,(r) 



a 



M(r 3 + e 3 M 3 )P 5 
P3 



=F 6aM(r 3 + e 3 M 3 ) 



P1+P2 
ft 



(4) 



with the functions Pi, ...Pq given in Appendix A of JP11. The upper and lower signs refer to 
co-rotating and counter-rotating particles, respectively. Compared to JP11, we modified Equ. 
to account for the facts that (Pi + P2), P3, and P5 can be negative and P3 should thus not be pulled 
out of the square root expressions. Over most of the allowed r-range, cr = 1; for some models P5 
has an extremum with P5 = at r a close to nsco an d = — 1 has to be chosen for r < r a to make 
L z smooth and differentiable. Equations Q and together with the equations 



Po 



and p 



-1 



(5) 



determine the four covariant (as well as the four contravariant) components of the momentum of 
the orbiting particles in GC. We restrict the discussion in this paper to prograde accretion disks. 

The simulations require to transform tangent vectors like the photon wave vector k from the 
GC system to the local inertial frame of the orbiting particles, called the Plasma Frame (PF) in the 
following. We define an orthonormal basis vector system (indices with a hat) to treat emission and 
scattering processes in the P F with one basis vector e^- chosen to be parallel to the four ye locity 
of the orbiting particles (e.g. Beckwith. Hawlev &i Krolik 2008 ; Shcherbakov Sz Huang! I20 111 ). The 
four basis vectors are given by the equations: 



dfi = G r / yj g rr , 

e = ae ( + ^ e f 

which define the transformation matrices e u „ through the relation 
a and (3 are the positive solutions of the equations e 4 - • e ? = and e 



(6) 



e; = 



The two constants 
1. We denote the 



inverse transformation matrices with a bar: e„ = < 
vector k then transform from the GC to the PF according to 



The components of a photon's wave 



(7) 



- 6 - 



and from the PF to the GC according to 

fc* = e" A ** (8) 



2.2. Radial Structure of the Thin Accretion Disk 

We consid e r sta ndard thin disk models with zero torque at the ISCO. For this case 



Page Sz Thornd ()1974l ) (called PT74 in the following) showed that mass, energy, and angular mo- 
mentum conservation alone determine the radial brightness profile of the accretion disk - exactly 
as in the case of a Shakura-Sunyaev accretion disk. PT74 derive the results for the general case of 
an axially symmetric metric given in the equatorial plane of the accretion disk in the form 

ds 2 = -e 2v dt 2 + e 2 ^(# - ujdt) 2 + e 2 ^dr 2 + dz 2 . (9) 

The conservations laws can then be used to show that the time average flux of radiant energy 
(energy per unit proper time and unit proper area) flowing out of the upper surface of the disk, as 
measured by an observer on the upper face who orbits with the time-average motion of the disk's 
matter, is given by the equation: 

F(r) = M° e - (u+ ^f(r) (10) 

with Mq being the radius independent time averaged rate at which rest mass flows inward through 
the disk. The function / depends on the momentum p of the orbiting particles: 

f(r) = ^r ^fdr (11) 
P<t> J nsco P 

where "," denotes ordinary partial differentiation and nsco is the radius of the ISCO. Comparison 
of Equation with Equation (pQ) for 9 = 90° and with z properly defined, gives the functions 
v, t/j, fj,, and u) as function of r. As p^r) and p^(r) can be inferred from Equations ([3][5]), it is 
straight forward to solve Equations (jlOp and (jlip numerically. We cross-checked the results by 
testing that the emitted luminosity integrated over the entire accretion disk equals Mq times the 
radiative efficiency (1 — -E(nsco))- 

The disk is assumed to have a temperature of 

( F \ 1/4 

T eS = — (12) 

V^SB/ 



with cjsb the Stefan Boltzman constant, and to emit a diluted blackbody spectrum with a hardening 
factor of /h = 1.8 (see below). 
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2.3. Photon Emission, Ray Tracing, and Scattering of the Polarized Photons 



Our code simulates the emission of photons from the accretion disk, the photon propagation 
through the spacetime, and the scattering of photons returning to the accretion disk. For each 
geodesic, the code keeps track of statistical information including for example the temperature of 
the disk segment which emitted the photon. This statistical information is used when the observed 
energy spectra are calculated. Although the discussion below uses the term photon, our treatment 
corresponds to the simulation of statistical ensembles of photon wave packages. In the following, 
x^ refers to the approximate location of a wave package, and k, EE, and f denote the mean wave 
vector, polarization degree and polarization v ector of the wave package, respectively (see MTW73, 
paragraph 22.5 and lGammie &: Leung] (|2012l ) for a discussion of the k and f). 



The numerical code simulates n photons emitted from the plane of the accretion disk (6 = 90°) 
in m radial bins logarithmically spaced from nsco to r max = 100, with n = 500 and m = 10, 000. 
Owing to the azimuthal symmetry of the problem, it is sufficient to simulate photons originating 
at (j) = 0. The photons are launched into the upper hemisphere (6 < 90°) with constant probability 
per solid angle in the PF, and with a n initial wave vector k p normalized such that Aj = |fe *| = 1 



in the PF. We use Table XXIV of (j Chandr asekharl Il960l . called C60 in the following) for the 
polarized emission from an optical thick atmosphere to calculate the statistical weight w cm for 
the chosen emission direction in the PF, and to calculate the initial polarization II of the photon. 
The polarization vector is initialized by setting /* = and choosing /* normalized to one and 
perpendicular to k % and (e^)\ After calculating the components of k and f in the PF, they are 
transformed into the GC system. The polarization degree is an invariant. 

The wave v ector k and polarizati o n vec tor f are parallel transported with a similar algorithm 
as described by iPsaltis Johannsenl (|2012l ) , extended to transport not only k but also f. The 
two Killing vectors £i = (1,0,0,0) and £2 = (0,0,0, 1) of the stationary axially symmetric metric 
imply the conservation of the photon energy and angular momentum at infinity: 



-kt 



-9tt 



dt 
dX 



9t<p 
d(j) dt 
^dX+^dX 



dX 



(13) 
(14) 



with the affine parameter A. These two equations are combined with dt/dX = fc* and d<ft/dX = k^ 
to calculate E~ and L 7 at the starting point of a geodesic (and also after the photon has been 
scattered). A fourth order Runge-Kutta algorithm is used to integrate the geodesic equation: 



d 2 x» 



dx a dx v 



dX' 2 ' av dX< dX< ' 
and to parallel transport the polarization vector according to the equation: 

Ax v 



dp 
dX' 



-p/i fa 
A auJ 



dX' 



(15) 



(16) 
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with Y^rjy denoting the Christoffel symbols and A' = E~A. 

Equations (|13M14p lead to the equations dt/dX' = {—g^ — bg t( p)/d and dcj)/dX' = (bgu + gt<f>)/d 
with b = Lj/E~ f and d = g^gtt — g t l which can be used to calculate dt/dX' and dcj)/dX' without the 
need for performin g the Runge-Kutta integ ration. Note that the complex-valued Walker-Penrose 
integral of motion ( Walker &: Penrose! 1970) tha t can be used to simplify the parallel transport of 



f through the Kerr metric ([Connors et al.l I198Q : IConnors Stark! [1980) was derived for vacuum 



solutions of Einstein's equation and cannot be used for the the metric of Equ. [TJ The accuracy of 
the integration is monitored with the invariants Ai = k 2 , A 2 = f 2 — 1, and A3 = k • f along the 
photon trajectory. 

We assume that photons can cross the equatorial plane between the event horizon and the IS CO 
without interacting with the matter plunging into the black hole. If a photon hits the accretion disk 
it is scattered using Chandrasekhar's formalism for scattering by a indefinitely deep atmosphere 
(C60, Section 70.3). After calculating the components of the wave and polarization vectors in the 
PF, the Stokes and Chandrasekhar parameters are computed. The direction of the scattered photon 
in the PF is drawn from a random distribution with equal probability per solid angle and Equation 
(164) and Table XXV from C60 are used to calculate the outgoing Chandrasekhar parameters. 
Referring to Fig. 8 of C60, a beam with Stokes parameters (Q = 0, U = 1) correspond to a 45° 
clockwise rotation of the polarization direction relative to the direction with (Q = 1, U = 0) when 
looking towards the origin of the coordinate system for both, the incoming and the outgoing beams. 

As an indefinitely deep atmosphere backscatters all incoming radiation and does not absorb 
any energy, it is clear that for each incoming photon direction given by the polar coordinates 
9q, (fio, the average statistical weight averaged over all outgoing photon directions 9, <fi should be 
unity. Chandrasekhar's formula gives the polarized components of the outgoing intensity / (energy 
emitted per unit time, accretion disk area, and solid angle) as function of the polarized components 
of the incoming energy flux ttF per unit area perpendicular to the direction of the incoming photon. 
The expression of the statistical weight for a certain scattering process reads: 

27T//J 2/iJ 

w sc = = — -. (17) 

HqttF fi F 

The factor of 2tt in the numerator of the second expression converts from "probability per solid 
angle" to "probability", the factor /i converts I into the outgoing energy flux per solid angle and 
accretion disk area, and the factor /io in the denominator converts ttF into the incoming energy flux 
per accretion disk area. We verified that averaging w sc over all directions in the upper hemisphere 
indeed gives unity. The simulation of a scattering event is completed by normalizing the temporal 
and spatial components of after scattering to the value before scattering (both in the PF), and 
calculating the GC components of k and f. 



- 9 - 



2.4. Analysis of the Simulated Events 

Photons are tracked until they reach a distance of r T = 10, 000 or until they get too close to the 
event horizon. Once a photon reaches r = r r , the components of the wave vector k r and the Stokes 
parameters are computed in the reference frame of a receiving observer with fixed coordinates. We 
present below most results in the coordinate frame of this observer with momentum p r = et/^/gu- 
The tangent vectors are given in terms of components (marked with a tilde) with regards to the 
orthonormal basis: 

e t = Pr = e t/y/9tt, (18) 
e^ = e r / \J g rr , 

and with e^ being a linear combination of e^ and e^ normalized to one and orthogonal to e^. 
Photons arriving in the lower hemisphere are mirrored into the upper hemisphere to compensate 
for the fact that the simulations cover only emission into the upper hemisphere. The process leaves 
the Stokes parameter Q invariant and flips the sign of U. 

Photons are accumulated in four ±4°-wide theta bins. Each photon contributes to the results 
with a statistical weight of: 

dN 

u> st = 27rAr— — — -w em W sc . (19) 
at dr d<p 

Here, Ar is the width of the radial bin from which a photon was launched, and dN/drdtdcj) is 
the number of photons emitted per GC unit time, radius, and azimuthal angle at the position 
of the radial bin. The factors w cm and W sc are the statistical weights for the emission process 
and all scattering processes, respectively. The latter equals the product of the weights w sc of all 
scattering events if the photon was scattered, and 1 i f not. We calculate d N/dr dt d(j) using a similar 



calculation as the one described in Appendix B2 of (jKulkarni et al.ll201ll ). The number of photons 

(20) 



emitted per area dA and time dt as measured in the PF is: 

dN F 



d?V <E> 

with d?V = dA dt. Assuming a diluted blackbody spectrum with a hardening factor of /h = 1.8, the 
mean PF energy per photon is <E> ~ 2.70 /h &b T c $. From the facts that the proper four-volume 
is an invariant and that the proper distance along the ^-direction is an invariant for boosts along 
the e^-direction, it follows that the proper 3- volume yj—gtr4> dt dr dcj) (with g tr( f, being the t-r-cj) part 
of the metric) is also an invariant. Taking into account that dN is an invariant and that d 3 V = 
^-gtr<j>dt dr d(j), we get: 

dN 



\/-gtr<t>—~— (21) 



dtdr dcp <E > 

with yj—gtr<f> = t for the Kerr metric and y/—gtr^ = (1 + M r )) r f° r the JP11 metric. 
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2-D maps of the emission are made using k®/k\ and kf/kl to give the photon direction in the 
image plane. Each simulated photon is used once, rotating a photon arriving at <f) to the azimuthal 
direction of the observer <p r = 0. No additional re-weighting is necessary for this approach, as 
one can see by considering that if a c/>-bin of a certain width A<p was adopted to select observed 
events, each simulated photon could be used N times with regularly spaced ^-offsets from to 
360° , choosing N such that the number of photons falling into the ^>-bin equalled unity. 

For all energy-resolved results, the photons contributes with a statistical weight of 

WE ^ = Wst J El ^J^i dE/ 1 ^~L dE (22) 

with eo = ft T e s k\ / k$ to a bin ranging from E\ to E2 ■ The factor k l r / k$ gives the redshift (or 
blueshift) from Doppler shifts following the emission of the photon and from the scattering(s) of the 
photon, and the gravitational redshift. We measure the polarization direction x (i- e - the direction 
of the electric field vector) from the projection of the spin axis of the black hole in the sky with 
X increasing for a clockwise rotation of the polarization direction when looking towards the black 
hole. 

The ray-tracing algorithm offers ample opportunity for consistency checks, i.e. k 2 = and 
f 2 = 1 can be checked along the photon trajectory. We compared the 2-D maps and energy spectra 
of the polarization degrees and polarization direction from our code with those from Figs. 1-3 o f 



(jSchnittman fc Krolikll2009l ) and found excellent agreement. Note that lSchnittman Sz Krolikl ((2009) 
adopt a different definition of the polarization vector. They use x = f° r a polarization parallel to 
the disk and x increases for a counter-clockwise direction of the polarization direction in the sky. 



3. Results 

We simulated the 10 parameter combinations listed in Table Q] (Models A- J). We assume a 10 
Mq black holes accreting at rates between 0.5 and 4 times 10 18 g s . The accretion rates were 
chosen to give a disk luminosity L-q = [1 — -E(nsco)] M of 10% of the Eddington luminosity. We 
considered Kerr black holes with spins of a = 0, 0.5, 0.9, and 0.99, and non-GR metrics with a = 0.5 
(e 3 = -30.6, -5, 2.5, and 6.3) and a = 0.99 (e 3 = -5, -2.5). Model E (a = 0.5, e 3 = -30.6) was 
chosen to have an ISCO nsco = 10 r g larger than that of a Schwarzschild black hole (nsco = 6r g ). 
Model H (a = 0.5, £3=6.3) was chosen to have exactly the same ISCO as the Kerr black hole with 
a = 0.5 (both nsco = 2.32 r g ). 

Figure Q] presents the event horizons in GC coordinates for all considered spacetimes. Note 
that the requirement of a closed horizon (and the absence of a naked singularity) limits the allowed 
range of £3 to £3 < e max with e max =6.75 for a = 0.5, and e max =0.01 for a = 0.99. For a = 0.99 and 
£3 = -5, the event horizon shows a kink close to the pole where two roots of the equation defining 
the event horizons cross. 
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Figures [2]|3] show the ray tracing images of the accretion disks for all simulated models for an 
inclination of 75° from the rotation axis. The lengths and orientations of the superimposed lines 
show the polarization degree and polarization direction, respectively. The brightness maps exhibit 
the well-known asymmetric appearance resulting from the relativistic beaming and de-beaming of 
the emission from the disk material approaching and receding from the observer with a velocity 
close to the speed of light, respectively, and the asymmetric gravitational lensing in the curved 
spacetime. The rings between the event horizon and the inner edge of the accretion disk come from 
photons orbiting around the black hole for multiples of ~ 180°. The brightness maps of the GR 
and non-GR models look somewhat similar. The largest differences result from a larger (smaller) 
ISCO of the non-GR models for negative (positive) 63-values. It is instructive to compare model 
A (e3 = 0, a = 0, Fig. [21 upper-left panel) with model I (63 = —5, o = 0.99, Fig. [U left panel): 
although the images of the ISCOs have similar angular diameters, model I shows the distortions 
that are typical for high spin values. From a theoretical viewpoint, the parameter M is the mass 
of the black hole as measured by a distant observer and is thus not a free parameter; the two 
images thus demonstrate that the metric from Equ. ([1]) with €3 7^ indeed leads to new observable 
characteristics, and does not merely correspond to the Kerr metric in different coordi nates. Note 



that in practice M may be known from the observations of the companion star, see e.g. (jOrosz et al 



20 111 ). For models with rather small ISCOs (Models D and H), the polarization direction of the 
emission from the inner disk exhibits a pattern that "circles around the black hole" owing to the 
dominance of the scattered emission. For the other models, the combination of GR effects and the 
competition between the direct and scattered emission results in a complex polarization pattern 
with subtle differences between the GR and the non-GR models. 

Figures [M3 present the accretion disk brightness F(r), the energy spectrum, the polarization 
degree spectrum, and the polarization direction spectrum for all simulated models, again for an 
inclination of 75° from the rotation axis. In the case of the Kerr black holes with fixed mass M 
(Fig. [5]), a determines the location of the event horizon and the ISCO. With increasing a, the ISCO 
moves towards smaller r, the disk brightness and disk temperature increase close to the black hole 
and the energy spectra extend to higher energies. The energy spectrum of the scattered radiation 
is harder than that of the direct emission owing to the facts that the scattered emission originates 
from the hot inner accretion disk and that the scattering can increase the photon energy owing 
to an additional Doppler boost. As a consequence of the hard spectrum and the relatively high 
polarization degree (compared to the direct emission), the scattered emission strongly impacts 
the polarization properties at higher energies. The polarization direction exhibits a swing from 
horizontal polarization at low energies (owing to the emission of the optically thick accretion disk) 
to vertical polarization at high energies (owing to scattered emission). The competition of the 
horizontal and vertical polarizations leads to a minimum of the polarization degree at intermediate 
energies. The swing of the polarization direction and the minimum of the polarization degree are 
observed at lower energies for the smaller ISCO models owing to the increased overall importance 
of scattered radiation. 
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The corresponding results are shown in Figs. [6] and [7] for the non-GR black holes with a=0.5 
and a=0.99, respectively. The results exemplify that the observational results strongly depend on 
the location of the ISCO: the smaller €3, the smaller risco> the harder the energy spectra, and the 
more pronounced is the energy dependence of the polarization properties. 

Under the assumptions made in this paper, GR and non-GR models produce qualitatively 
different energy spectra and polarization spectra for 63-values which lead to nsco > 6r g . As an 
example, Model E (with a = 0.5, £3 = —30.6, risco = 10 r g ) exhibits a softer energy spectrum 
and less variation of the polarization degree and polarization direction than any of the Kerr-models. 
Existing X-ray spectroscopic data can already be used to exclude such negative 63-values. 

The difference between GR and non-GR models is rather small for models with the same nsco- 
As an example, Fig. [8] shows GR and non-GR models with approximately the same risccr va hies. 
The GR and non-GR models produce almost indistinguishable energy spectra and similar - but 
not identical - polarization properties. The differences between the GR and non-GR models are 
somewhat larger for rapidly spinning black holes because more photons are emitted and propagate 
close to the event horizon where the JP11 metric deviates most strongly from the Kerr metric. 
Although the GR model C and the non-GR model H have the same ISCO and show almost identical 
flux energy spectra, they exhibit different polarization energy spectra. The latter depend more 
strongly on the underlying spacetime as the competition of the direct and scattered emission leads 
to more pronounced observational signatures owing to the very different polarization properties of 
the two emission components. 

4. Summary and Discussion 

In this paper we explore the possibility to test GR in the strong gravity regime with X- 
ray spectropolarimetric observations of black holes in X-ray binaries. In the thermal state, the 
accretions disks, and the emission and scattering geometry and processes are relatively simple and 
well understood, making these systems attractive for the test of the underlying spacetime. We 
developed a code to simulate the polarized emission of the accretion disk. The code computes the 
radial structure and radial brightness profile of the accretion disk for GR and non-GR spacetimes. 
Furthermore, it parallel transports the wave vector and the polarization vector and accounts for 
scattering based on the classical results of Chandrasekhar (C60). 

We used the code to study the observational differences between Kerr black holes and the black 
holes described by the phenomenological metric of JP11. The main effect of the JP11 metric for 
X-ray spectral and polarization measurements is to allow for other ISCOs than those predicted by 
GR for a certain M and a combination. The X-ray spectropolarimetric observations will allow us 
to measure allowed and forbidden regions in the plane of the parameters a (black hole spin) and 
€3 (quantifying the deviation from the Kerr metric). However, for large regions of the parameter 
space, the approximate degeneracy between the black hole spin a and the deviation parameter €3 
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will make it difficult to distinguish between GR and non-GR models. 



The GEMS mission achieves the best sensitivity (~1% polarization degree f or a mCrab source ) 
with a photoelectric effect polarimeter operating over the 2-10 keV energy range (jBlack et al.ll2010i ). 
GEMS would be able to detect the polarization of the X-ray emission from bright X-ray binaries 
like Cyg X-l and GRS 1915+105 with an excellent signal to noise ratio. For a deep (~ 10 6 sec) ob- 
servation of a source with a flux of 300 mCrab, statistical and systematics errors of the polarization 
in 1 keV wide bins would be of the order of 0.2% (Kallman et al., private communication). The 
errors would be sufficien tly small to allow us to dist inguish between the Kerr models shown in Fig. 
[5] (see the discussion by ISchnittman Krolild |2010| . for more details). However, the GEMS sensi- 
tivity would not be sufficient to distinguish between the GR and non-GR models shown in Fig. [HJ 
In practice, systematic errors on the model predictions associated with remaining uncertainties of 
the structure of the accretion disk and its atmosphere (influencing e. g. the extent of de polarization 
due to Faraday rotation, and the spectral hardening factor /h, see iDavis et al J 12009 ). as well as 
systematics associated with the subtraction of "contaminating" emission components would further 
complicate this task. Extending the polarimetric coverage to lower and higher energies would help 
to further constrain models. Polarimetry de t ectors with a broader bandpas s have been discussed 



in the literature fe.g. IWeisskopf et al. 



2010 



Krawczvnski et al 



2011 



2006 



Lei et al. 



1997; 



Bellazini et al. 



2010 



Marshall et al 



2012J). Excellent sensitivity below 2 keV would make it possible to 



constrain the impac t of Faraday rotat ion and thus the magnetic field structure in the accretion disk 
and its atmosphere ([Davis et al.ll2009l ) . A polarimetry mission with excellent sensitivity in the hard 
X-ray band (e.g. the BEST mission, see iKrawczynski et all l2012l ) would allow us to characterize 
the polarization properties of the harder emission components and to reduce the errors associated 
with subtracting these components. Other metrics may imply different observational effects, and 
it is clear that it would be desirable to explore the space of non-GR metrics in a more systematic 
fashion. 

Testing fundamental physics with astronomical observations requires a proper understanding 
of the astrophysics of the observed systems. In our case, we need to continue to refine accretion 
disk and radiation models before studies like the ones presented in this paper can be used to 
test GR. As mentioned above, numerical models are now used to test the more than forty year 
old conjecture of a vanishing torque at the ISCO, and give more detailed information about the 
physical properties of accretion disks (e.g. about th e radius at which the Compton optical depth 
becomes smaller than unity, see iKrolik et all 120051 ). Even though numerical models have made 
impressive progress, they are still missing important physics including detailed modeling of the effect 
of radiative energy transport. Future studies of the predicted polarization signatures could improve 
on the simplified modeling of the emission and scattering processes adopted in this paper. Such 
calculations could employ a more detailed model of the accretion disk with information about the 
structure of the plasma and the magnetic field. More accurate modeling of the emission, scattering, 
and absorption processes could account for the composition and ionization state of t he disk material , 



for different atomic transitions, multiple scattering, and Compton scattering (see iNagirner 



1962 
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Loskutov fc Sobolevlll98ll. 1198a Davis Hubenvl 12006) . Recently, iBroderick k Blandfordl (|200l 
2004 ). and lGammie &: Leung (|2012l ) developed methods for modeling the propagation of polarized 
rays through magnetized plasmas in curved spacetimes which may be used to model the polarization 
of the emission from turbulent and magnetized accretion disks accounting for the Faraday rotation 
of the polarization direction. 



Johannsen &: Psaltid (|2012l ) studied the possibility to use Fe K-a fluorescent line observations 



of black holes in X-ray binaries to distinguish between the Kerr metric and the metric of Equ. ([I]). 
Similarly as in this paper, they find that the parameters a and 63 are degenerate for some regions 
of the parameter space. For these regions, distinguishing between the different metrics will require 
observations with a very high signal to noise ratio, an exquisite understanding of the accretion disk 
properties, and excellent control over systematics associated with the subtraction of the underlying 
continuum emission. 



Bambil (|2012l ) discusses the possibility to use the observed jet power to get an independent 
handle on t he black hole spin. Assum ing that black hole jets are powered by the Blandford-Znajek 
mechanism ([Blandford Znajekill977l ). the correlation between the jet power and the black hole 
spin can be used to break the degeneracy between a and €3. Although this technique has promise, 
its practical utility will be limited in the foreseeable future by our imperfect understanding of the 
jet launching mechanism and systematic effects, e.g. the correlation between the magnetic field 
close to the black hole event horizon and the black hole spin. 

The temporal variabilit y of the obser v ed flux es and polarization properties offer additional 
diagnostics. For example, Horak fc Karas ( 2Q06al Jbh study observational signatures produced by 
moving clouds that scatter emission and emphasize that gravitational lensing can result in time de- 
lays for parts of the signal. It may be possible to use time resolved spectropolarimetric observations 
of such transient events to constrain the geometry of the system and the underlaying spacetime. 
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Table 1. Parameters describing the simulated metrics 



Model 


M[M ] 


MflO^gs- 1 ] 


a 


£3 


r H [r g ] 


nsco [r g ] 


A 


10 


2.45 








2 


6 


B 


10 


1.7 


0.5 





1.87 


4.23 


C 


10 


0.90 


0.9 





1.44 


2.32 


D 


10 


0.53 


0.99 





1.14 


1.45 


E 


10 


4.00 


0.5 


-30.61 


3.08-3.13 


10.00 


F 


10 


2.33 


0.5 


-5 


1.87-1.96 


5.79 


G 


10 


1.27 


0.5 


2.5 


1.80-1.87 


3.28 


H 


10 


0.88 


0.5 


6.33 


1.74-1.87 


2.32 


I 


10 


1.88 


0.99 


-5 


1.22-1.87 


5.09 


J 


10 


1.49 


0.99 


-2.5 


1.14-1.71 


2.61 
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Fig. 1. — The outer event horizons of the considered black holes. The upper left panel shows Kerr 
black holes with (from the outside inwards) a = 0, 0.5 , 0.9, and 0.99; the upper right panel shows 
black holes with a = 0.5 and (from the outside inwards) 63 = —30.6, —5,0,2.5 and 6.3; the lower 
panel shows black holes with a = 0.99 and (from the outside inwards) €3 = —5, —2.5, and 0. 



-21 - 




-W [r /r r ] -kf/kj [r /r r ] 




Fig. 2. — Simulated images (0.1-10 keV) of four accreting Kerr black holes for an inclination of 
75°. The results are shown for a = (top-left, Model A), a = 0.5 (top-right, Model B), a = 0.9 
(bottom-left, Model C), and a = 0.99 (bottom-right, Model D). See Tabled] for all the parameters 
of the different models. The logarithmic color scale shows the 0.1-10 keV photon number flux. The 
length of the bars show the polarization degree and the orientation shows the polarization direction 
of the electric field vector. Far away from the black hole the polarization degree is ~4%. 
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Fig. 3.— Same as Fig. O but for a = 0.5 and e 3 = -30.6 (top-left, Model E), a = 0.5 and e 3 = -5 
(top-right, Model F), a = 0.5 and e 3 = 2.5 (bottom-left, Model G), and a = 0.5 and e 3 = 6.3 
(bottom-right, Model H). 
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Fig. 4. — Same as Fig. O but for a = 0.99 and e% = —5 (left panel, Model I) and a = 0.99 and 
£3 = —2.5 (right panel, Model J). 
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Fig. 5. — Accretion disk brightness in the plasma frame (top-left), energy spectrum (top-right), 
polarization degree (bottom-left) and polarization direction (bottom-right) for four Kerr black 
holes with a = (solid black line, Model A), a = 0.5 (dashed blue line, Model B), a = 0.9 (dotted 
green line, Model C), and a = 0.99 (dash-dotted red line, Model D). The inclination is 75°. The 
polarization direction x 1S measured from the projection of the spin axis of the black hole in the 
sky with x increasing for a clockwise rotation of the polarization direction when looking into the 
beam. 
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Fig. 6. — The same as Fig. [5] but for black holes with a = 0.5 and €3 = —30.6 (dotted black line, 
Model E), a = 0.5 and €3 = —5 (dash-dotted magenta line, Model F), a = 0.5 and £3 = (dashed 
blue line, Model B), a = 0.5 and £3 = 2.5 (solid green line, Model G), and a = 0.5 and €3 = 6.3 
(dashed red line, Model H). 
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Fig. 7. — The same as Fig. [5]but for black holes with a = 0.99 and £3 = —5 (solid black line, Model 
I), a = 0.99 and €3 = —2.5 (dashed blue line, Model J), and a = 0.99 and £3 = (dash-dotted red 
line, Model D). 
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Fig. 8. — The same as Fig. [5] but for black holes with a = and €3 = (solid black line, Model A), 
a = 0.5 and €3 = —5 (dash-dotted magenta line, Model F), a = 0.5 and 63 = 6.3 (dashed red line, 
Model H), and a = 0.9 and €3 = (dotted green line, Model C). 



